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ABSTRACT: We study the evolution of a cosmological baryon asymmetry produced via leptogen- 
esis by means of the full classical Boltzmann equations, without the assumption of kinetic equi- 
librium and including all quantum statistical factors. Beginning with the full mode equations, we 
derive the usual equations of motion for the right-handed neutrino number density and integrated 
lepton asymmetry, and show explicitly the impact of each assumption on these quantities. For the 
first time, we investigate also the effects of scattering of the right-handed neutrino with the top 
quark to leading order in the Yukawa couplings by means of the full Boltzmann equations. We find 
that in our full Boltzmann treatment the final lepton asymmetry can be suppressed by as much as a 
factor of ~ 1.5 in the weak wash-out regime (K < 1), compared to the usual integrated approach 
which assumes kinetic equilibrium and neglects quantum statistics. This suppression is in contrast 
with the enhancement seen in some previous studies that considered only decay and inverse decay 
of the right-handed neutrino. However, this suppression quickly decreases as we increase K. In the 
strong wash-out regime (K > 1), the full Boltzmann treatment and the integrated approach give 
nearly identical final lepton asymmetries (within 10% of each other at K > 3). Finally, we show 
that the opposing effects of quantum statistics on decays/inverse decays and the scattering processes 
tend to reduce the net importance of scattering on leptogenesis in the full treatment compared to 
the integrated approach. 
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1. Introduction 

Leptogenesis [1] provides an attractive explanation for the baryon asymmetry of the universe. The 
ingredients of this mechanism are simple: heavy right-handed Majorana neutrinos are added to the 
standard model of particle physics, which on the one hand explains naturally the smallness of the 
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observed neutrino masses by way of the see-saw mechanism [2]. On the other hand, the out-of- 
equilibrium decay of the heavy neutrino states into leptons and Higgs particles violates CP, from 
whence a lepton asymmetry can be generated. This lepton asymmetry is then partially transformed 
into a baryon asymmetry by anomalous processes of the standard model called sphalerons [3]. 

In the recent past a huge step forward has been made towards understanding in detail the 
processes of leptogenesis. Relevant studies include leptogenesis in a supersymmetric context [4], 
thermal effects [5, 6], analytic formulae for the final efficiency factor [7], the role of flavour [8, 
9, 10, 11], as well as leptogenesis with CP violation coming only from the measurable low- 
scale Pontecorvo-Maki-Nakagawa-Sakata matrix [12]. Furthermore, it has been pointed out that 
the classical Boltzmann equations are insufficient to describe the transition region between the 
flavoured and the unflavoured regimes [13]; a full quantum-mechanical description in terms of 
density matrices is necessary. On this front, the quantum-mechanical Kadanoff-Baym equations 
have been investigated for toy models in extreme out-of-equilibrium situations [14, 15]. 

On a different front, the classical Boltzmann equations have been solved for the first time 
for single momentum modes [16]. As one of Sakharov's conditions [17], departure from thermal 
equilibrium is crucial for the dynamic creation of a baryon asymmetry. In the leptogenesis sce- 
nario, out-of-equilibrium conditions are achieved when interactions are no longer able to maintain 
the momentum distribution function of the right-handed neutrino at its equilibrium value as the 
universe expands. To simplify the calculation, this non-equilibrium process is traditionally studied 
by means of the integrated Boltzmann equations [6, 18, 19, 20], whereby the equations of motion 
for the distribution functions of all particle species involved are integrated over momentum such 
that only the evolution of the number densities is tracked. However, in order for the integrated 
equations to be in a closed form, it is necessary to neglect quantum statistical behaviours (e.g., 
Pauli blocking) and assume kinetic equilibrium for all particle species, including the right-handed 
neutrino. For particle species with gauge interactions these assumptions seem justifiable. For the 
right-handed neutrino however, their validity is not immediately obvious. 

To estimate the effects of kinetic equilibrium and quantum statistics, the Boltzmann equations 
for the individual momentum modes have been solved in [16], taking into account only the decay 
and inverse decay of the right-handed neutrino within the unflavoured framework. More recently, 
the mode equations have been used to study the effect of a pre-existing asymmetry and the soft 
leptogenesis scenario, again including only decays and inverse decays [21]. In the present work, 
we extend on these previous studies by considering also scattering processes of the right-handed 
neutrino with the top quark. 

The paper is organised as follows: We describe the basic set-up of our scenario in section 2. In 
section 3 we study the simplest scenario involving only the decay and inverse decay of the right- 
handed neutrino. In section 4 we include scattering processes mediated by the Yukawa interaction 
of heavy neutrinos with the top quark. We conclude in section 5. 

2. Basic set-up 

We concentrate on the simplest case of "vanilla-leptogenesis", in which a lepton asymmetry is es- 
tablished from the decay and scattering of the lightest heavy right-handed neutrino N\ . We neglect 
the decay of the two heavier neutrino states [22], assuming that any lepton asymmetry pro- 
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duced from these decays will be efficiently washed out by the N\ interactions. Therefore we will 
drop the subscript "1", and refer to the lightest right-handed neutrino simply as N in the following. 
Furthermore we will work in the one-flavour approximation, since flavour effects do not change 
the kinetic consideration for the mode equations. 

Generally, the Boltzmann equation (BE) for a right-handed neutrino (RHN) in a Friedman- 
Lemaitre-Robertson-Walker framework can be written as 

- \ PN \H ^ = C D [f N ] + C s [f N ] , (2.1) 
at o\Vn\ 

where /jv is the phase space distribution of the RHN, p jy the RHN momentum, and H the Hubble 
parameter. On the right-hand side, the collision integrals Cd [/jv] and Cs [/jv] encode respectively 
the interactions of the RHN due to decays into leptons and Higgs (D) and scattering processes via 
Yukawa interactions with the top quark (S). 

The Boltzmann equation for leptons (anti-leptons) with phase space distribution fi (fj) has a 
similar form to Eq. (2.1), save for the replacements //y — > fi (/jy — > fj) and pjv ^ P/ (Pjv ^ Pj)- 
Since we are interested in the asymmetry between leptons and antileptons, it is convenient to define 

fl-l = fi~ h (2-2) 
and the corresponding Boltzmann equation 

% 1_|P ' |F §^ = ° D [fl ^ + ° S [/ '"' ] ' (2 " 3) 
where C D ,s [fi-j] = C D ,s [fi] ~ C DtS [fj] ■ Integrating /,_ T over the lepton phase space, i.e., 



n 



= d3pi fl ~~ v {2A) 



with gi = 2, gives us the lepton asymmetry per co-moving photon, 



Nt_j = (2.5) 



where n^- = (((3)/7r 2 )g 1 T 3 , with g 7 = 2, is the equilibrium photon density. 

Equations (2.1) and (2.3) can be further simplified by transforming to the dimensionless co- 
ordinates z = M/T and y« = |p«|/T, where M is the mass of the RHN [23]. Using the relation 
dT/dt = —HT, the differential operator dt — \pi\Hd\ p .\ becomes zHd z , and consequently 

()]s( '- !,) ' (C D [f N (z,y)]+C s [fN(z,y)}), (2.6) 



and 



dz H(M) 

dfi--i(z,y) _ z 
dz H(M) 



{CDif^y^+Csif^y)}) , (2.7) 



with H(M) = y/4:ir 3 g*/45(M/Mpi), where M P1 = 1.221 x 10 19 GeV is the Planck mass, and 
g* = 106.75 corresponds to the number of relativistic degrees of freedom in the standard model at 
temperatures higher than the electroweak scale. 
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The Boltzmann equations (2.6) and (2.7) encode how a lepton asymmetry is generated and 
washed out in an expanding universe given some specific particle interactions. Here, it is useful to 
define an efficiency factor that measures the amount of asymmetry that has survived the competitive 
production and wash-out processes, 

/c = | e" 1 JV,.!- (2-8) 

The quantity e quantifies the amount of CP violation in the interactions, while N t _j is the lepton 
asymmetry produced. In the limit of a vanishing wash-out and a thermal initial abundance for the 
RHN, the efficiency factor has a final value Kf = 1. 

3. Decay and inverse decay 

We consider first the simplest possible scenario of thermal leptogenesis, in which only the decay 
and inverse decay of the RHN into lepton I and Higgs $ pairs contribute to the evolution of /at, 
i.e., we set Cs = in Eqs. (2.6) and (2.7). The decay and inverse decay of the RHN violate CP 
through interference of the tree-level and the one-loop diagrams. 

The collision integral for the RHN in the decay-inverse decay picture has the following form: 

Cd \f N ] =^r- / J-^r, (2vr) 4 5 4 ( PN -pi- p») 



2E N J 2^(2vr) 3 2£*(27r) 3 
x [h fi (1 - In) (\M^n\ 2 + \M^ N \ 2 ) (3-1) 
-f N (1 - /,) (1 + U) {\M N ^i\ 2 + \M N ^j\ 2 ) ] , 

where Ei and pi are, respectively, the energy and 4-momenta of the particle species i, and Ma 
denotes the matrix element for the process A. At tree level, the squared matrix element summed 
over all internal degrees of freedom for the decay of the RHN into a pair of lepton and Higgs 
particles is given by 

\M N ^ = 2 {m ° m 2 D)ll p lPN , (3.2) 

where v = 174 GeV is the vacuum expectation value of the Higgs particle, and the Dirac mass 
matrix mo is connected to the Yukawa coupling matrix X u via mo = v\ v . 
The integral (3.1) can be readily reduced to a one dimensional form [24] 

Mr rf r(E N +VN)/2 

ENPN J{E n - Pn )/2 

where 

rhiM 2 



Cd [In] = -pf, / dp* - f N ) - Ml - fi)(l + /*)] , 

ENPN J(Em-vm)/2 



(3.3) 



Frf = ^ (3 " 4) 
is the total decay rate in the RHN's rest frame, with 

mi = (mi C )n (3.5) 
M 

the effective neutrino mass [20], to be compared with the equilibrium neutrino mass 

£ — ,2 



16^2^/7/ v 2 

m * = 7= T7 — • (3-6) 

3^5 Mpi 
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The decay parameter 

K = mW)~^ ( } 

controls whether the RHN decays in equilibrium (K > 1) or out of equilibrium (K < 1). 

For leptons participating in the same decay and inverse decay processes, the collisional integral 
is given by 

n r,, 1 /" d 3 p N d 3 p^ 4 4 

x [/iv (1-/0 (1 + /$) \M N ^*i\ 2 (3 - 8) 

- (1-/at) \M^ n \ 2 ]. 

An analogous expression for the anti-leptons can be derived by replacing fi — > /j, jv^$; — > 
•Mjy^j, and M^i^n -^^j^n- Some useful relations exist between the matrix elements 
following from CPT-invariance [18]: 



\M N ^\ 2 = \M^ N \ 2 = |M | 2 (1 + e), 
\M N ^i\ 2 = \M^n\ 2 = \M \ 2 (1 - e), 



where | .Mo 1 2 is the tree level matrix element given in Eq. (3.2). 

The collision integral (3.8) suffers from the problem that a lepton asymmetry is produced even 
in thermal equilibrium. This can be remedied by including contributions from the resonant part of 
the AL = 2 scattering process l& <-> 7$ [6, 18]. We implement this remedy following the method 
developed in [7], and add to the collision integral (3.8) the term 

Ufj (1 - f N )\M^ N \ 2 uh - hh (1 - f N )\M^ N \ 2 uh , (3.10) 

where 

l-A^ivlsub = \M\ 2 AL=2 - e \M \ 2 , 
\M^ N \ 2 uh = \M\ 2 AL=2 + e \M \ 2 , 

where \M\ 2 AL=2 is negligible for M <C 10 14 GeV [25]. 1 

In the following subsections, we review first the derivation of the conventional integrated 
Boltzmann equations, which neglects quantum statistics and assumes kinetic equilibrium for the 
RHN. We then remove step by step these assumptions, in order to examine their effects on the 
efficiency factor re. The scenarios to be examined and their associated assumptions are summarised 
in table 1 . 

3.1 Case Dl: Integrated Boltzmann equations 

In the integrated approach conventionally used in the literature [6, 7, 18], the time evolution of 
number densities are tracked in favour of the phase space distributions /j. This is achieved 
by integrating the Boltzmann equations (2.6) and (2.7) over momentum. However, the integrated 
equations have no closed forms unless we make certain simplifying assumptions: First, we ne- 
glect factors stemming from Pauli blocking for fermions and induced emission for bosons, i.e., 



'Reference [21] includes terms in addition to Eq. (3.10) in order to avoid asymmetry production in thermal equilib- 
rium. However, the same analysis also shows that the quantitative difference between this and our approach is negligible. 
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Table 1: Scenarios considered in the decay/inverse decay picture and their associated assumptions. Case Dl 
corresponds to the conventional integrated Boltzmann approach, while Case D4 was previously investigated 
by Basb0ll and Hannestad [16]. 





Assumption of kinetic equilibrium 


Including quantum statistics 


Section 


CaseDl 


Yes 


No 


3.1 


Case D2 


No 


No 


3.2 


Case D3 


Yes 


Yes 


3.3 


Case D4 


No 


Yes 


3.4 



we approximate 1 ± fi ps 1 [18]. Second, all standard model particles are taken to be in ther- 
mal equilibrium due to their gauge interactions and their distribution functions approximated by a 
Maxwell-Boltzmann distribution, f^ q = e~ E ^ T . 

With these assumptions and using energy conservation, we find 

Ml = e -^ +El )/T = e -E N /T = ; eq ^ 

so that the collision integral (3.3) simplifies to 

C D [f N ] = — / dp* [ffi - f N ] . (3.13) 

&NPN J(E N - PN )/2 

Integrating (3.13) over p$ and inserting into Eq. (2.6), the BE for the RHN distribution function 
becomes 

df N zF ri M ( . 

-dT~ H(M)E N V»- fN >- (3 " 14) 
To make further inroads, we assume kinetic equilibrium holds for the RHN, i.e., its distribution 
function /at can be expressed as $n / f°N ~ un/u^, where is the RHN number density. Then 
one can easily integrate Eq. (3.14) over the RHN phase space to obtain 

where K = T r{ /H(M) (cf. Eq. (3.7)), and r if (M/E N ) = {F r{ /n eq ) f d 3 PN /{2ir) 3 Q (M/E N ) 
is the thermal average of the decay rate [18]. The thermally averaged dilation factor is given by 
the ratio of the modified Bessel functions of the second kind of first and second order, (M/En) = 
K 1 {z)/K 2 (z). 

Dividing Eq. (3.15) by the equilibrium photon density n^ q , we obtain the Boltzmann equation 
for the quantity Nn = riN/n^ q [7], 

°*L = -D(N N -N$), (3.16) 

with 

D -«*(£). <M7> 

and 

N eq (z) = -z 2 K 2 (z). (3.18) 
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Here, an inconsistency in the integrated approach is visible: all particles, i.e., N, I and $, are 
assumed to follow the Maxwell-Boltzmann distribution function. However, when calculating N 1 ^ 1 , 
we must use a Fermi-Dirac distribution for the RHN, = [3 ((3) g N T M 2 /(8vr 2 )] K 2 (z), with 
C(3) « 1.202 and gjy = 2, in order to reproduce a realistic equilibrium RHN to photon density 
ratio. This leads to an extra prefactor (3/4)£(3) in the definition of N 1 ^ compared to a strictly 
Maxwell-Boltzmann approach. 

For the lepton asymmetry, Eqs. (2.7), (3.8) and (3.10) combine to give the BE for the lepton 
distribution functions, 



dfl- 



dz 



z 2 K 



-iyf 



(3.19) 



where Em = Em IT. Using energy conservation and assuming kinetic equilibrium for the RHN, 
Eq. (3.19) can be integrated over um to give [16] 



dfl-l 



z 2 K - z +4! " 
■ e 



UN 



n 



N 



n 



cq 
N 



dz 2yf 

We further assume that both kinetic chemical equilibrium prevails for the leptons such that 

/£. = e -(^~M)/T _ e -{El+n)/T „ 2 (/u/T) e - El /T t 

fl-l 



(3.20) 



(3.21) 



n 



n 



cq 



with chemical potential \i <C 1 and n^ q the lepton equilibrium number density. Thus, integrating 
over the lepton phase space, we obtain the equation of motion for the number density 



dn, 



z 3 KT 3 



dz 



2tt' 2 



KUz) 



n 



l-l 



n 



cq 



2e 



ec i\ i 
nN -nM 



n 



cq 
N 



(3.22) 



where K\{z) is the modified Bessel function of first kind. Following [7] we rewrite Eq. (3.22) in 
terms of the lepton asymmetry per co-moving photon , 



dN, 



-i 



dz 



-W ID N l _ i + eD(N N -NP), 



where 



1 1 N cq 

W ID ^-Kz 3 K l (z) = -D J ^ 



(3.23) 



(3.24) 



quantifies the strength of the wash-out due to inverse decays, and A^ eq = 3/4. Note that, as with 
the RHN, when evaluating iV, eq it is necessary to use a Fermi-Dirac distribution for the leptons, 
n^ q = (3/4) (C(3)/7r 2 )^T 3 , with gi = 2, to ensure a realistic lepton to photon density ratio. 

Figure 1 shows the final efficiency factor Kf, defined in Eq. (2.8), as a function of K for several 
different initial RHN abundances [26]. 
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Figure 1: Final efficiency factor for different scenarios of thermal and non-thermal leptogenesis. Shown are 
k f for a thermal (dashed/red), a vanishing (dot-dot-dash/light green), and several cases of dominant initial 
RHN abundance. A dominant initial abundance is realized if a scalar field responsible for inflation decays 
exclusively into the RHN, which then dominates the energy density of the universe. The coupling strength 
between the scalar field and the RHN can be connected to an energy scale: 10 9 GeV (solid/blue), 5 x 10 s GeV 
(dashed/dark green), 3.75 x 10 s GeV (dot-dash/purple) and 10 8 GeV (dot-dash-dash/taupe). 



3.2 Case D2: Dropping the assumption of kinetic equilibrium 

Since the RHN is very heavy — its mass scale corresponds to the temperature of the thermal bath 
during the period of leptogenesis — it is not a priori obvious that decays and inverse decays would 
occur fast enough to establish kinetic equilibrium. Thus the assumption of kinetic equilibrium 
for the RHN might lead to sizable deviations from an exact treatment. In this section we drop this 
assumption in our calculation of the efficiency factor. We retain however our other assumptions: that 
all equilibrium distribution functions are of the Maxwell-Boltzmann form, and quantum statistical 
factors are negligible. 

Dropping the assumption of kinetic equilibrium for the RHN means that it is now necessary to 
solve Eq. (3.14), rewritten here as 

^ = (e~ £N ~ In) , (3-25) 

OZ ON 

individually for all possible values of the dimensionless RHN energy £ pf. For the calculation of the 
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lepton asymmetry, the relevant equation is Eq. (3.19) which we reproduce here 

- ' 



dz 



DC 

— iy 
4 «f 



2 dy N ^ [Uf l _ 1 -2e(f N -f e I ?)] 



(3.26) 



Again, this equation must be solved for all possible values of the lepton momentum yi, and the 
resulting /,_j(y;) summed according to Eq. (2.5) to give N Alternatively, using energy con- 
servation and assuming kinetic and chemical equilibrium for the standard model particles, we can 
integrate Eq. (3.26) over the lepton phase space to obtain a single equation of motion for iV ; _j, 



dN, 



l-l 



z 2 K 



dz 



POO P 

L dyi h 



, v dVN^g [N l _- l Q-2e{f N -Q)]. 



(3.27) 



We find the second approach to yield more stable results. 

3.3 Case D3: Boltzmann equations with quantum statistical factors 

In Case D3 we reinstate Pauli blocking factors for fermions and factors due to induced emission for 
bosons, but adopt again the assumption of kinetic equilibrium for the RHN. Consistency requires 
that we use the Fermi-Dirac and the Bose-Einstein distribution functions respectively for fermions 
and bosons in thermal equilibrium, instead of the classical Maxwell-Boltzmann distribution func- 
tion. 

With these assumptions in mind, we integrate the collision integral (3.3) over pq, to obtain the 
BE for the RHN, 



Of. 



N 



z 2 K njy 



n 



cq 
N 



1 



dz £ N y N 



n 



cq 
N 



v£n + l 



log 



sinh ((£ N - y N ) /2) 
sinh {(£ N + y N ) /2) 



(3.28) 



where we have used /at//^ 5 = (1 + e £N )fN ~ n^/n c ^. Integrating over the RHN phase space 
and normalising to the photon number density yields 



dN 



N 



K 



(Nn-N c «) f 
Jo 



dy 



dz K 2 {z) J "» N S N + 1 

We note that the integral over the RHN phase space has no simple analytic form. Therefore it 
remains necessary to perform the integration numerically. 

The BE for the lepton asymmetry including all quantum statistical factors and assuming kinetic 
equilibrium for all particle species has the following form: 



VN 



1 



log 



sinh((£jv -y N )/2) 
smh((£ N +y N )/2) 



(3.29) 



dz 



z 2 K r°° 

~2yJ Jfz^l 

4,y t 



2 dy N TT- 



£n 



(h + ^ /£)(/,_! + eF+) -2e% Q{1 + f 9 ) 



n 



N 



n 



N 



(3.30) 

where F + = fi + fj ~ 2 / ; eq . After integrating over the lepton phase space and normalising to the 
photon number density we arrive at 



dN, 



i-i 



dz 



z 2 K 
4 

CO 



f 

Jo 



dyi 



-iy 



, VN 
2 dy N — 



AfCq 
N 



N 



3 



2e 



fP 



2£ -TTeq fj^ (1 

N 



(3.31) 



with Nf? given in Eq. (3.18). 
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3.4 Case D4: Complete mode equations 

Here we include all statistical factors and make no assumption of kinetic equilibrium for the RHN. 
Integrating Eq. (3.3) over p$ gives the BE for the RHN, 



df N 

dz 



K 



£nun 



-l + f N + e £ »f N 
e £jv + i 



log 



sinh((gjy -y N )/2) 
sinh((£Ar + yjv)/2). 



(3.32) 



The equation for the lepton asymmetry in this case is similar to Eq. (3.30), except we do not assume 
kinetic equilibrium for the RHN, 



dfi-i 

dz 



= "^r Z 2_ iv? dy N !g [(/* + f N M-i + eF+) - 2ef N (1 + /*; 



(3.33) 



Integrating over the RHN momentum yields 



dz 



1 



4y { 



/•oo 

■ / dyi 
Jo 



a y N 



(U + In) (-N l _- l + 2e) fP-2ef N (1 + /*) 



(3.34) 



where we have assumed, as usual, thermal equilibrium for the standard model particles. 

3.5 Results and discussions 
3.5.1 Right-handed neutrino 

Figure 2 shows the time evolution of the comoving number densities of the RHN for the four 
different cases described above, assuming a vanishing initial RHN abundance. We have picked two 
values for the decay parameter: (i) K = 0.1, lying in the weak wash-out regime, is shown on the 
left hand side, and (ii) K = 10, lying in the strong wash-out regime, is shown on the right hand 
side. For reference, we also plot the time evolution of the RHN equilibrium number density Nff. 

The general behaviour of the RHN abundance evolution is similar for all four cases. In the 
weak wash-out regime, there is a net production of RHN by inverse decays at high temperatures 
z < 1. At z ~ 4, the RHN abundance overshoots the equilibrium density and continues to grow 
until z ~ 5, when a net destruction of RHN by decays into l& pairs begins to push its abundance 
slowly back down to the equilibrium value. Equilibrium is reached finally at z ~ 20, beyond 
which the RHN abundance falls off exponentially with z, as expected for all non-relativistic particle 
species in thermal equilibrium. Contrastingly, in the strong wash-out regime, the stronger coupling 
brings the RHN abundance to its equilibrium value already at z ~ 1. Its subsequent evolution is 
then simply governed by equilibrium statistics: at z ~ 4 the RHN becomes nonrelativistic and 
hence its abundance is suppressed by exp(— M/T). 

In both the weak and strong wash-out regimes, the difference between Cases Dl and D2, 
which exclude quantum statistical factors, and their counterparts Cases D3 and D4, which include 
quantum statistics, is most visible at z < 1. The RHN abundance is almost an order of magnitude 
larger in the latter two cases than in the former. This is because during the high temperature RHN 
production phase, using the correct Bose-Einstein equilibrium distribution function for the Higgs 



-10- 



K = 0.1 



K = 10 




Figure 2: Time evolution of the comoving RHN number density Nm and of the absolute value of the lep- 
ton asymmetry N t _j, assuming two different coupling strengths K — 0.1, 10, and e — 10~ 6 . The four 
scenarios within the decay-inverse decay only framework are shown: Solid/red line denotes Case Dl, long 
dashed/green D2, dotted/blue D3, and short dashed/magenta D4. See table 1 for a short summary of each 
scenario. For reference, we also indicate the equilibrium RHN number density N 1 ^ 1 in dot-dash/cyan. 

boson /$ substantially enlarges the phase space available for the inverse decay process &l — > N at 
low E$. This effect is far stronger than the phase space restriction due to Pauli blocking by the final 
state RHN, as can be seen from the phase space factors in Eq. (3.3). As the temperature drops and 
the RHN becomes nonrelativistic, the effects of quantum statistics also diminish, since kinematics 
now prevents the low energy <I> and I states from contributing to the collision integrals. 

Interestingly, the assumption of kinetic equilibrium leads to no visible effects in either the 
weak or strong wash-out regime. Comparing Cases Dl and D2 (both assume Maxwell-Boltzmann 
statistics), their RHN abundances are virtually identical. The same is true for Cases D3 and D4, 
which include quantum statistical factors. 
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0.01 0.1 1 10 100 

K 

Figure 3: The final efficiency factor k as a function of K for the four scenarios within the decay-inverse de- 
cay picture, assuming a vanishing initial RHN abundance. Solid/red line denotes Case Dl, long dashed/green 
D2, dotted/blue D3, and short dashed/magenta D4. 

3.5.2 Lepton asymmetry 

The time evolution of the corresponding absolute value of the lepton asymmetry is shown in the 
lower panel of Figure 2. A negative lepton asymmetry is produced at high temperatures by RHN 
production from inverse decays. At z ~ 5 in the weak wash-out regime (z ~ 1 if strong wash-out), 
decays come to dominate over inverse decays, thus reversing the direction of the asymmetry pro- 
duction, and eventually flipping the sign of the asymmetry to positive. When the RHN abundance 
begins to fall off exponentially, the asymmetry also asymptotes to a final, constant value. 

In the weak wash-out regime (K = 0.1) the asymmetries produced in Cases D3 and D4 
which include quantum statistical factors are always larger in magnitude than those produced in 
Cases Dl and D2 which assume Maxwell-Boltzmann statistics throughout the whole temperature 
range considered. The change of sign also occurs slightly earlier in D3 and D4. These effects 
can be understood as follows. From, e.g., Eq. (3.33), we see that the production of a negative 
lepton asymmetry at high temperatures by inverse RHN decays is significantly enhanced when 
we take proper account of the Bose-Einstein statistics for the Higgs boson. Like the case of the 
RHN abundance, this effect dominates over the phase space suppression due to the Fermi-Dirac 
statistics of the lepton and the RHN. As we progress to lower temperatures, RHN decays begin to 
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dominate over inverse decays, thereby reversing the direction of the lepton asymmetry evolution. 
Since quantum statistics speeds up RHN production and brings its abundance up to the equilibrium 
threshold earlier, the transition from decay to inverse decay domination — and hence the turning 
point in the asymmetry evolution — also happens earlier. As a result, the asymmetry flips sign a 
little earlier in Cases D3 and D4 than in Dl and D2, and has more time to grow to a larger positive 
value before the exponential fall-off of the RHN abundance shuts down the asymmetry production. 

In the strong wash-out regime (K = 10), a similar behaviour is also visible at z < 1. As 
we progress to lower temperatures, however, Cases D3 and D4 end up producing less asymmetry 
than Cases Dl and D2. This is because for K > 1, the wash-out rate plays a dominant role in 
determining the final asymmetry. Here, quantum statistics enlarges the phase space of the wash-out 
term from in Eq. (3.26) to (/$ + fN)fi_j in Eq. (3.33), thus forcing the lepton asymmetry 

to flip sign even earlier, and continuing on to dampen it to a slightly smaller positive value. 

Again, as with the RHN, the assumption of kinetic equilibrium has virtually no effect on the 
asymmetry evolution: the differences between Cases Dl and D2, and between Cases D3 and D4 
are generally at the percent level, too small to be visible in Figure 2. 

Finally, Figure 3 summarises the lepton asymmetry produced in the four cases, in terms of the 
final efficiency factor Kf defined in Eq. (2.8), as a function of the decay parameter K. For all values 
of K considered, the assumption of kinetic equilibrium can be seen to produce a minute (< 5%) 
difference in Kf between Cases Dl and D2 and between Cases D3 and D4. Quantum statistics, on 
the other hand, has a generally stronger effect on the final lepton asymmetry. In the weak wash-out 
regime (K < 1), inclusion of quantum statistical factors (Cases D3 and D4) enhances Kf by a 
factor of ~ 1.5 relative to Cases Dl and D2 which assume Maxwell-Boltzmann statistics. In the 
strong wash-out regime (K > 1), the effect of quantum statistics is to suppress Kf by up to 20% at 
K ~ 10, but reduces to the percent level at K ~ 100. 

4. Scattering processes 

In this section we enlarge our picture of thermal leptogenesis to include tree level scattering pro- 
cesses of the RHN with the top quark , e.g., Nl — ► qt, which are of O (/?|A 2 ). These interactions 
lead to an additional production channel for the RHN and contribute to the wash-out processes. 
Until now these scattering processes have only been considered using the integrated Boltzmann 
equations [6, 7]. Here we provide for the first time a solution of the full set of Boltzmann equations 
at the mode level. 

We do not consider interactions with gauge bosons, nor include CP violation in 2 — > 2 or 
1(2) — ► 3 processes, which are of higher order in the Yukawa couplings. CP violation from these 
processes was considered in [9, 27, 28, 29], where it was shown that at high temperatures CP 
violation from scattering is the main source of lepton asymmetry production. However the final 
asymmetry depends also on the strength of the wash-out processes; It turns out that in the weak 
wash-out regime (K < 1) CP violation in the scattering processes tends to suppress the asymmetry 
production, while in the transition (K ~ 1) and strong (K > 1) wash-out regimes its contribution 
is small to negligible. 
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Table 2: Scenarios including scattering with the top quark and their associated assumptions. Case SI cor- 
responds to the conventional integrated Boltzmann approach, while Case S2 involves solving the full set of 
Boltzmann equations at the mode level. 





Assumption of kinetic equilibrium 


Including quantum statistics 


Section 


Case SI 


Yes 


No 


4.1 


Case S2 


No 


Yes 


4.2 



In the following we first recall the treatment of scattering processes in the integrated picture, 
before we proceed to write down the full set of mode equations including the relevant scattering 
terms. Table 2 summarises the assumptions of these two scenarios. 

4.1 Case SI: Scattering in the integrated picture 

Analogous to Eqs. (3.16) and (3.23), the integrated Boltzmann equations including scattering take 
the following form [7], 

^ =-{D + S)(N N -N$), (4.1) 

9N l l / eaN 

_ti =- W N l _ 1 + eD (N N - N$) , (4.2) 

where S accounts for the production of RHN from scattering processes, and the wash-out rate W 
contains also a contribution from these processes. The scattering rate S itself consists of two terms, 
S = 2 S s + 4 St, coming respectively from scattering in the s-channel and in the i-channel. One 
factor of 2 stems from contribution from processes involving anti-particles, and another factor of 2 
in the t-channel term originates from the u-channel diagram. 
In general, the scattering rates are defined as 



Q — ^ s ' t 

H z 

where 



M 1. 



s,t 



s.t 



(4.3) 



(4.4) 



24((3)g N 7T 2 K 2 (z)zf 

Note that an additional factor of 4/ (3) C(3) appears in this definition compared to the definition of 
reference [7]. This is due to the Fermi-Dirac statistics used in our derivation. The quantity l Sjt is 
an integral 

I, tt = j J (W a a , t (9) V^Ki (y/^j (4.5) 
of the reduced cross-section a St t, given by [4] 



3/i t 2 AT mi , . 
<?s,t = -r 1 - — Xs,t (x) , (4.6) 

47T V Z 

where x = ^/z 2 , and h t = fh(T) is the top Yukawa coupling, to be evaluated at the relevant 
energy scale (or temperature) T by solving the renormalisation group equation. See Appendix B. 



-14- 



The functions Xs,t{%) are defined as 



V,(.r) (^) , (4-7) 



Xt{x) = X 



x 



x -2 + 2a h 1 - 2dh , { x -1 + a h 
H r- lo g 



(4.8) 



where we have introduced ah = m<$,/M as an infrared cut-off for the t-channel diagram, and m$ 
is the mass of the Higgs boson which presumably receives contributions from interactions with the 
thermal bath, i.e., its value does not correspond to that potentially measured at the LHC. The value 
of m$ can in principle be deduced from a thermal field theoretic treatment of leptogenesis, and the 
analysis of [6] found m$(T) ~ 0.4 T. However some open questions still remain and hence in the 
present work we prefer to adopt a value of = 1CT 5 , used first by Luty in [30]. 

It is convenient to rewrite the s- and £-channel scattering rates S Sj t in terms of the functions 
f S: t defined as 

J~«x«,t Of/; 2 ) v-pki (W) 

'•■'^ = m^) — — • (4 - 9) 

such that 

Ss,t = ^jfs,t, (4.10) 

and the total scattering rate is given by 

S=§^(fs(z) + 2f t (z)), (4.11) 



where 



with rhi given by Eqs. (3.5), and [7] 



K s = % (4.12) 

mi 



where is defined in in Eq. (3.6). 

Since the scattering processes with the top quark change the lepton number by one unit, they 
contribute also to the wash-out of the asymmetry. The total wash-out rate is given by 

W = W ID + W AL=1 , (4.14) 

where Wid denotes the contribution from inverse decay defined in Eq. (3.24), and Wal=i from 
scattering in the s- and t- channels, 

W A L=i = W s + 2W t , (4.15) 

with 

w _N N ri _N% N N 
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and 



N' 



Hz N° q 



cq 

N St. 



(4.17) 



The lepton scattering rates are given by T l s t = N^/N^ r Sjt . Using Eq. (3.24), the two contribu- 
tions Wid and Wal=i are related by 



W AL=1 = 2W ID 



( N N 



15 [l^ Ss + 2St] ' 



so that 



W = W ID 



1 + d{ 2 ^ 3 ^^ 



(4.18) 



(4.19) 



gives the total wash-out rate. 

Figure 4 shows the various rates D, S, W and Wid as functions of z assuming K = 0.1. For 
other choices of K, the corresponding rates evolve with z in a similar fashion, but with magnitudes 
scaling with K. 
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4.2 Case S2: Complete mode equations including scattering 

The basic BE for the distribution function of the RHN is given by 

H(M) df N 



dz 



C D [In] + 2 Cs, s [f N ] + 4 C s ,t [In] , (4.20) 



where again one factor of 2 stems from contribution from processes involving anti-particles, and 
another factor of 2 in the i-channel term originates from the u-channel diagram. The decay-inverse 
decay collision integral Cd is given in Eq. (3.1), the s-channel scattering integral is 



i=l,q,t 

X [(1 - f N )(l - fOftfg - f N fi(l - f t )(l - f q )\ , 



(4-21) 



and a similar expression exists for the i-channel scattering integral Cs,t [In], but with the appro- 
priate matrix element Mt, and the replacements fi <-> f q . 
The analogous equation for the lepton asymmetry is 

^ = C D [/ M ] + 2 Cs, s [ fl _j\ + 4 C s , t [U , (4.22) 
where Cd [fi-j] = Cd [fi] — Cd [fj] can be constructed from Eq. (3.8), and 

csAfi-i]=^f n 7j^.(^ 4 (pi+PN- Pq -Pt)\M S \ 2 

2hlJ i=N, q A > 1 ( 4 - 23 ) 

(fN (ft + fq-l)~ ftfq) ■ 

Replacing M s with M. t and f q <-> /jy in Eq. (4.23) yields the integral Cs,t [fi-j] ■ 

The collision integrals (4.21) and (4.23) are nine-dimensional and can be reduced analytically 
down to two dimensions. We give in appendix A the final reduced integrals and, as an example, 
the full reduction procedure applied to the s-channel collisional integral for tracking the RHN 
following the method of [31, 32]. A general treatment of scattering kernels in kinetic equations can 
be found in [33]. 

4.3 Results and discussions 

4.3.1 Scattering vs decay-inverse decay 

Figure 5 shows snapshots of the RHN distribution function in Case S2 relative to an equilibrium 
Fermi-Dirac distribution at time z = 0.2, 1, 5, as well as the RHN number density normalised to 
its equilibrium value as a function of z for an interaction strength of K = 0.1. These are compared 
with their counterparts assuming decay and inverse decay only (Case D4). Figure 6 is similar, 
except for K = 10. Clearly, including scattering processes speeds up the equilibration of the RHN 
distribution function, especially at high temperatures (z < 1). This effect is more significant for 
small values of K, since for large K values decays and inverse decays are already fast enough to 
establish equilibrium. 
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Figure 5: Snapshots of the RHN distribution function Jn/^n at 2 — 0-2, 1, 5, and the RHN abundance 
Nn/N^ as a function of z, assuming K = 0.1. Solid/red line denotes Case Dl, long dashed/green D4, 
dotted/blue SI, and short dashed/magenta S2. See tables 1 and 2 for a summary of the scenarios. 



Looking at the time evolution of the RHN number density we see a corresponding increase in 
Npf at high temperatures when scattering is included (Case S2), compared to the case with decays 
and inverse decays only (Case D4). The equilibrium density is also reached at an earlier time (or 
higher temperature). The integrated approach shows a similar behaviour, with Case S 1 predicting 
a large RHN abundance at high temperatures and hence faster equilibration than Case Dl. 

Figure 7 shows the time evolution of the lepton asymmetry, again for the two characteristic 
values of the decay parameter K = 0.1, 10. As discussed earlier, we have explicitly ignored CP 
violation in the scattering processes, so that they have no direct influence on the lepton asymmetry. 
This assumption is the reason why, in the weak wash-out regime (K = 0.1), the asymmetry evolu- 
tion at high temperatures (z < 1) in Case S2 is virtually identical to that in its decay-inverse decay 
only counterpart Case D4. The same behaviour can also be seen when comparing Cases SI and 
Dl. Here, decays and inverse decays of the RHN alone source the creation of a lepton asymmetry. 
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Figure 6: Same as Figure 5, but for K = 10. 



Since for K < 1 the asymmetry evolution at high temperatures hinges primarily on inverse decays 
and is as yet unaffected by such external factors as the RHN abundance and wash-out processes, 
the inclusion of scattering processes has no visible effect on N,-,. 

However, scattering can still affect the asymmetry production in two indirect and competing 
ways: (i) the larger RHN abundance produced via scattering processes at high temperatures forces 
the lepton asymmetry to flip sign earlier, thereby generating a larger positive lepton asymmetry, 
and (ii) scattering leads to additional wash-out of the lepton asymmetry The first effect dominates 
for coupling strengths lying in the weak wash-out regime (K < 1), eventually leading to a larger 
asymmetry in Cases SI and S2, compared with their decay-inverse decay only counterparts Dl 
and D4 as shown in Figure 7. For stronger couplings (K > 1), the second effect dominates; in fact, 
Figure 7 shows that the additional wash-out due to scattering suppresses the lepton asymmetry 
production in Cases SI and S2 already at high temperatures z < 1, compared with the decay- 
inverse decay only scenarios D 1 and D4. 
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Figure 7: Time evolution of the absolute value of the lepton asymmetry | JV, j | for two different cou- 
pling strengths K. Shown are the two cases including scattering processes SI (dotted/blue) and S2 (short 
dashed/magenta), and two scenarios Dl (solid/red) and D4 (long dashed/green) within the decays-inverse 
decay only framework. For reference we also plot the RHN equilibrium abundance (dot-dash/cyan). 



4.3.2 Complete treatment vs integrated approach 

The complete treatment differs from the integrated approach in that in the latter case we assume ki- 
netic equilibrium for the RHN and neglect all quantum statistical factors. As we saw in section 3.5, 
the assumption of kinetic equilibrium tends to underestimate by a tiny amount the RHN abundance 
at z < 1. This can be understood from Figures 5 and 6 as a result of the more efficient production 
of low momentum RHN states, which in turn contribute more to the momentum integral. 

Quantum statistics, on the other hand, has very different effects on the scattering and the 
decay-inverse decay collision terms. As we saw in section 3.5, in the decay-inverse decay sce- 
nario, quantum statistics always enhances the interaction rates through the enlarged Higgs boson 
phase space density /$ at low E$. For the scattering processes, since all participants are fermions, 
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the role of quantum statistics is to reduce the phase space and hence suppress the interaction rates. 
In general, however, we expect quantum statistics to be more important for decay/inverse decay 
than for scattering. This is because in the decay-inverse decay case the enhanced phase space due 
to /$ at low E$ can in principle be infinite, while Pauli blocking for fermions participating in 
scattering, e.g., 1 — //, suppresses the phase space by at most a factor of 1/2. 

The difference between the RHN abundance and the lepton asymmetry evolution in the com- 
plete and the integrated treatments can then be understood in terms of a competition between the 
three aforementioned effects. 

Consider first the RHN abundance. At z <C 1 the dominant RHN production channel is scat- 
tering. Here, suppression of the production rate due to quantum statistical factors competes with 
the small enhancement due to our dropping the assumption of kinetic equilibrium. The net result 
is that both Cases SI and S2 give very similar RHN abundances as shown in Figures 5 to 7. At 
z ~ 0.3, decay/inverse decay becomes comparable to scattering (see Figure 4). Here, the enhanced 
decay rate due to quantum statistics in Case S2 pushes up RHN production relative to Case SI. This 
effect is more prominent in the weak wash-out regime than in the strong wash-out region, since in 
the former case the RHN abundance is further away from equilibrium. Progressing further in z, we 
see that the RHN abundances in Cases SI and S2 become virtually identical already before reach- 
ing the equilibrium value. This is in stark contrast with the decay-inverse decay only scenarios, 
where the RHN abundances in Cases Dl and D4 clearly cross the equilibrium threshold at different 
times. 

Consider now the evolution of the lepton asymmetry (lower panel of Figure 7). Comparing 
Cases SI and S2 in the weak wash-out regime (K = 0.1), quantum statistics in the latter scenario 
enhances the production of a negative lepton asymmetry at high temperatures. This effect is due 
solely to phase space enhancements in the inverse decay term, since we have assumed explicitly that 
scattering does not violate CP. At z ~ 4, the production of lepton asymmetry reverses direction 
as RHN decays begin to dominate over inverse decays. As mentioned earlier, quantum statistics 
causes this reversal to happen earlier in the decay-inverse decay only scenario by bringing the RHN 
abundance to the equilibrium threshold at an earlier time. When including scattering, however, the 
RHN abundances in both Cases SI and S2 cross the equilibrium threshold at almost the same time, 
as discussed in the previous paragraph. This means that the evolution of their corresponding lepton 
asymmetries also turns around at roughly the same time. Since at the time of the turn-around Case 
S2 has a more negative asymmetry than Case SI, the net effect is that the asymmetry in Case S2 
flips sign at a slightly later time than in Case SI, and subsequently grows to a smaller positive 
value. 

The effects of quantum statistics on the lepton asymmetry evolution in the strong wash-out 
regime (K = 10) can be similarly understood, except that we must consider also the role of the 
wash-out terms. At z < 1, the wash-out rate is dominated by scattering. However, as shown in 
Figure 4, decay/inverse decay becomes comparable to scattering at z ~ 0.3 and is the dominant 
wash-out process at z > 1. Thus, from z ~ 1 onwards, the net effect of quantum statistics is 
to enhance the wash-out rate. This effect can be seen at the turn-around of the lepton asymmetry 
evolution: the stronger wash-out rate in Case S2 forces the lepton asymmetry evolution to reverse 
direction at a slightly earlier time than in Case SI. However, since at the time of the turn-around 
Case S2 has a more negative asymmetry than SI, the asymmetries in both cases end up flipping 
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Figure 8: The final efficiency factor k / without and with scattering terms: Dl (solid), D4 (long dashed), SI 
(dotted), and S2 (short dashed). 

signs at almost the same time and grow to nearly identical values. 

Finally, Figure 8 shows the final efficiency factors as a function of K for the integrated ap- 
proach and the complete mode treatment, both including and excluding scattering. For Cases SI 
and S2 which include scattering, we note that their difference is rather large in the weak wash- 
out regime (K < 1), with the integrated approach overestimating kj by up to a factor ~ 1.5 at 
K ~ 0.01 compared to solving the complete mode equations. But this difference decreases as we 
increase K. At K > 3 the integrated approach underestimates Kf by less than ~ 10%. 

It is also interesting to note that the relative contribution of scattering processes to the final 
efficiency factor is smaller in the complete mode calculation than in the integrated approach. In 
the weak wash-out regime, including scattering enhances the final efficiency factor from decays 
and inverse decays by up to a factor of ~ 30 in the integrated scenario. In the complete mode 
calculation, however, the enhancement is only a factor of ~ 15. Similarly, in the strong wash-out 
regime, scattering reduces kj by up to 20% in the integrated picture, compared to below 10% in 
the complete treatment. 



5. Conclusions 

In this paper we have for the first time studied leptogenesis by means of the full Boltzmann equa- 
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tions incorporating all quantum statistical terms without the assumption of kinetic equilibrium, and 
including scatterings of the right-handed neutrino with quarks. This is of particular relevance for 
the creation of the cosmological baryon asymmetry due to the required deviation from thermal 
equilibrium and the energy dependence of all interactions. As the simplest possible set-up to study 
these effects, we have considered only an asymmetry being created by the lightest right-handed 
neutrino and have neglected potentially important flavour effects. 

In the conventional approach, i.e., neglecting quantum statistical factors and assuming kinetic 
equilibrium, considering only decays and inverse decays is known to give a rather precise approx- 
imation of the final baryon asymmetry in the interesting strong wash-out regime. Hence, we start 
by considering this case, which also offers the possibility to study the influence of various effects 
separately and in detail. 

Interestingly enough, dropping the assumption of kinetic equilibrium has almost no effect 
on the evolution of the RHN number density and the lepton asymmetry. Taking the full energy 
dependence of interactions into account changes the final efficiency factor by 5% at the most. 

Including all quantum statistical factors has somewhat larger effects. These factors tend to 
enlarge the phase space available for neutrino production by inverse decays, thus significantly 
boosting the RHN abundance and the "wrong-sign" asymmetry being produced at high temper- 
atures. Further, they lead to an earlier domination of decays over inverse decays, thus speeding up 
the production of the final asymmetry. In the weak wash-out regime this leaves more time for the 
production of an asymmetry, thus leading to a boost in the final asymmetry by ~50%. 

In the strong wash-out regime, on the other hand, the final asymmetry is suppressed by between 
20% at K = 10 and 1% at K = 100. This is again due to the enlarged phase space of inverse decay 
processes which act as wash-out terms, thus reducing the asymmetry compared to the case where 
quantum statistical factors are neglected. 

In the case of scatterings of RHNs off quarks, in contrast to decays and inverse decays, quan- 
tum statistical factors reduce the phase space available, since all external particles in these processes 
are fermions. Hence, quantum effects generally tend to reduce the importance of these scatterings. 

Nonetheless, at high temperatures (z < 1) scattering processes increase the amount of RHNs 
being produced, thus making leptogenesis more efficient. On the other hand, at low temperatures 
(z > 1) they act as wash-out terms thereby reducing the produced asymmetry. The first effect, 
i.e., the more efficient production of RHNs dominates in the weak wash-out regime, thus leading 
to a larger final lepton asymmetry compared to the case where only decays and inverse decays 
are included. In the strong wash-out regime the effect from increased wash-out dominates, i.e., 
including scatterings leads to a somewhat reduced asymmetry here. 

This is qualitatively in line with results obtained in the integrated picture, i.e., neglecting quan- 
tum statistical factors and assuming kinetic equilibrium. However, since quantum statistical factor 
enhance decays and inverse decays while suppressing scatterings of RHNs with quarks, the net 
influence of these scattering processes is reduced when quantum factors are included. This signifi- 
cantly reduces the spread of results for the final efficiency factor, particularly in the weak wash-out 
regime. 



-23- 



Acknowledgments 

We thank Steen Hannestad, Josef Pradler, Max Huber and Georg Raffelt for useful discussions. 



Appendices 

A. Reduction of the scattering collision integrals 
A.l s-channel 

A.l.l Right-handed neutrino 



The full collision term for the s-channel in Eq. (4.20) is 

_ d 3 pi 
(1 

i=l,q,t 



CsMn] = ^ I II j^^(^Hpn +Pi-Pt- Pq )\M s \ 2 A^HfNjt, f t ,f q ), 



with phase space factor given by 



(A.l) 



Af } (fN,fi,ft,f q ) = [(1 - /iv)(l - fi)ftf q - f N fi(l - ft)(l - /,)] • (A.2) 

The matrix element M s is summed over all internal degrees of freedom of the particles in the initial 
and final states, including colour and isospin, and is given by 

|^| 2 = 24/ l ?^™^i, (A.3) 

where fh\ = (jn^mnj /M is the effective neutrino mass [20], v = 174 GeV the vacuum expec- 
tation value of the Higgs field and h 2 the top Yukawa coupling given in Appendix B. 
We work in the centre-of-mass frame, i.e., 

PN + Pi = Pt + Pq = q- (A.4) 

In general, the 4- vector delta function can dealt with using the relation 



S (p? - Mf) = S (El - (| Pi | 2 + M^) 



+ 



(A.5) 



v 2^| Pi | 2 + M 2 2^| Pj | 2 + M 2 ) 

Using this relation, and the fact that we consider all particles except the RHN to be massless, i.e., 

El,t,q = |PM,q|; 



E n = V\pn\ 2 + M 2 , 



(A.6) 
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we can integrate over the quark energy, 

d 3 p, 



I ^ 64{PN + Pi-Pt-P^ = J dE « d " P« M \J PgD ©(^) S(E N + E l -E t - E q ) 



@(E N + Ei- E t ) 



x 5 3 (ptv + Pi - Pt - p,) 
= 5 (E N + Ei-E t - \pn + Pi- Pt\) 

2 I PAT + Pi ~ Pt\ 

= HEn + e , -b .-| q -p.|) e(EN + El _ Et) 

2|q-Pt| 

= 5 ((E N + Ei- E t ) 2 - |q - p t | 2 ) Q(E N + Ei- E t ). 

(A.7) 

Similarly, we can rewrite 

= J dEi 6 (Ef - \pi\ 2 ) 0(£,) d 3 Pl (A.8) 
= J dE l 6(Ef-\q-p N \ 2 ) e(£,)A, 

where the last equality follows from changing the variable from pi to q = pjy + P;, and hence d 3 pi 
to d 3 q, and the integration is over the lepton energy E\. 
We choose an explicit coordinate system, 

q=|q| (0,0,1), 

Pat = | PN I (0, sin ry, cos rf) , (A.9) 
p t = E t (cos(f) sin i?, sin ^> sin cost?), 

and obtain the following quantities: 

s = (p* + Plf = (Pt+p q ) 2 = (E N + Ei) 2 - |q| 2 , 
s-M 2 



PNPl = Y 

PqPt = |, 

|2 _ |„|2 , I |2 „ „ I 1 2 , r 2 



(A. 10) 



|q " Pt| = Iql + IPtl - 2q ■ pt = |q| 2 + £?t - 2|q| £ t cos tf, 
|q - Piv| 2 = |q| 2 + IpaH 2 - 2q • pat = |q| 2 + |ptv| 2 - 2 |q| |pjv| cost/. 



The matrix element in these coordinates reads 



\M S \ 2 = 6h 2 ^ (E N + El ) 2 -M 2 -H 
v 2 {Ex + Etf-H* 
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and the delta functions are given by 

6 ((E N + Ei- E t f - |q - p t | 2 ) = S ((E N +Ei- E t f - |q| 2 - E 2 t + 2|q| E t cost?) 

- -J— 5 (cos 7? - E't-jEN + E.-E^ + H^ 

^(^-Iq-pjvl 2 ) = < 5(^ 2 -|q| 2 -|p iv | 2 + 2|q||p J v| cos 77) 

1 . ( _ | Pj y| 2 -^ 2 + |qn 
2|q||Piv| V OS1] 2\ q \\ PN \ )' 

(A. 12) 

Collecting all terms we get for the collision integral 

C sAfN] = 0TP \ , 5 / ^dcostidtfi^dEtdE^q-r- -i -_L— |.M S | 2 

2E N (2iry> J 4vr 2E t 2|q| |paH 2|q| ' 

/ | Pj y| 2 -^ 2 + |q|^ ./ , E 2 -(^ + J B i v-^) 2 + |q| : 
x d cos 77 — : cos v - 



2|q||Piv| 7 V " 2|q|£7, 
x Af ) (/at, /,, / t , /,) 0(£j) G(E t ) Q(Ei + E N - E t ) 

^ d cos t? d cos r) dE t dE x d\q\ |.M S | 2 



W^IPivl 



,. | PJV | 2 -£ 2 + |qH J £;2_ (E +Ejv _ jE;) 2 + | j: 

X COS 77 j— j-p j COS 77 — 



2|q||Piv| ) \ 2HE t j 

x Af ) (f N , /,, f t , f q ) e(E,) e(E t ) e(s, + ^ - 

(A. 13) 

Here, in the first equality, we take d 3 pt = E% dE t d cos ■& d(f), and average over the direction of 
the incoming RHN by integrating over dfi7v/(47r), where dfJjv = dO d cos 77, because of rotational 
invariance (cf. Eq (A.9)). In the second equality, we take d 3 q = Air |q| 2 <i|q|, and integrate over all 
azimuthal angles. 

The two remaining angles 7? and 77 run in the range 

cos 7?, cos 77 G [—1,1]. (A. 14) 

Since apart from the delta functions the integrand does not depend on either angle, integrating over 
these ranges effectively lead to new integration limits for the (/-integral: 

cost?=1 => q€[E l + E N ,-E l -E N + 2E t ], 
cos 7? = -1 => q G [-£7, - £jv, E t + E N - 2E t ], 
cos?7=l =>• ^ G [--B/ +PN,E t +Pn], 
cos?7 = -l ^ € [—Ei — pjv, ^ — Pn\i 

where = |pjv|, and q = |q|. Putting these conditions together we get 

sup [\2E t - E x - E N \, \Ei - p N \] < q < inf [E t + E N ,E L +p N ] . (A. 16) 

Since Ei + E^ > Ei + p^ this reduces to 

sup [\2E t -Ei- E N \, \E t - p N \] ^q^Ei+pN. (A.17) 



(A. 15) 
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Thus, the integration over cos 77 and cos •& effectively gives rise to a combination of 9 functions in 
the remaining 3-dimensional integral. Together with existing G functions in Eq. (A. 13), we define 

Q = e(q-\2Et-Ei-E N \) Q(q-\E l -p N \) 6 (£7, + p N - q) Q(E t + E N - E t ) (A.18) 

to collectively denote all 6 functions appearing in the remaining integral. Note that we have omitted 
writing out @(Ei)@(E t ), since positive particle energies are understood. 
Next, we use the relations 

6 (q - \2E t -E x - E N \) = 1-9 (\2E t - E x - E N \ - q) , (A.19) 

and 

9 (\2E t -Ei- E N \ - q) 9 {Ei + p N - q) = 9 (\2E t -E x - E N \ - q) , (A.20) 

the latter following from the fact that in Eq. (A.18) vanishes unless \2E t — E\ — E^\ < E\ +pn- 
With these we split the function ft into two parts (i.e., U = Qi + ^2): 

^1 = 6(^ + ^-^)6(^-1^-^1) @(Ei+ PN -q), (A.21) 
n 2 = - &(Ei + E N - E t ) e(q- \E t - p N \) 6 (\2E t - E l - E N \ - q) . (A.22) 

Equation (A.21) can be further split into two parts at 6 (q — \Ei — pn\) using the relation 

@(Ei- PN ) + e(p N -Ei) = l, (A.23) 

from which we find 

fila = e{E t + E N - E t ) e(q- (Et - p N )) 8 (Et + p N - q) 6 (E l - p N ) , (A.24) 

ni 6 = e(Et + e n - E t ) e( q - ( PN - s,)) e (E t + PN - q ) e ( PN - Et) , (a.25) 

so that Oi = Qi a + fin,. 

Similarly, Eq. (A.22) can be split at 9 (q — \Ei — pn\) into two parts, $^2 = ^2a + ^2b, by 
way of the relation (A.23): 

n 2a = ~ Q(Ei + E N - E t ) Q(q- (Et - p N )) 8 (\2E t - E l - E N \ - q) 9 (E t - p N ) , 

(A.26) 

n 2b = - e(E t +E N - E t ) Q(q- (p N - E t )) 9 (\2E t -Ei-E N \-q) Q(p N - E t ) . 

(A.27) 

One further split is possible at 9 (\2E t — Ei — E^\ — q) using 

9 (Ei + E N - 2E t ) + 9 (2E t -E x - E N ) = 1. (A.28) 
Putting this relation in Eqs. (A.26) and (A.27), we find the combinations 

9 (q - (Ei - p N )) Q(Ei + E N -2E t -q)^ E t <^ (E N + p N ) , 

Q(q- (Ei - p N )) Q(2E t -Ei-E N -q)^ E t >\ (2E t + (E N - p N )) , 

{ (A.29) 
© (q ~ (PN ~ Ei)) 0(Ei + E N -2E t -q)^ E t <- (2E t + (E N - p N )) , 

e(q- (p N - Ei)) Q(2E t -E N -Ei-q)^E t >^ (E N + p N ) , 
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with which we can write down the four parts of n 2 : 

Sha,i = -@(^(E N + p N ) - E^j G (q - (Ei - PN )) (Ei + E N - 2E t - q) 

x@(Ei-p N ), (A.30) 
^2a,u = ~ @(Ei + E N - Et) 6 (E t - ^(2Ei + E N - p N )^j 9 (2E t - E x - E N - q) 

x&(q-(E l - p N )) Q(E t - p N ), (A.31) 
^2b,i = - 6 Q(2^ + E N - p N ) - E^j 9 (Ei + E N - 2E t - q) 

X G (q - (p N - E^) Q(p N - Ei) , (A.32) 
n 2 b,u = ~ Q(Ei + E N - E t ) 6 (E t -\(E N + PN )^j 9 (2E t - E x - E N - q) 

x@(q-(p N - Ei)) Q(p N - Ei), (A.33) 

such that Q 2a = ^2a,i + ^2a,u and Q 2 b = ^2b,i + ^2b,u- 
Finally, collecting all terms we obtain the relation 

^ = ^2 = + + + ^2a,ii + ^2b,i + ^2b,ii, (A.34) 

so that the remaining 3-dimensional integration in the collision integral (A. 13) can be equivalently 
written as 

CsAfN] = Y\ o7 3 I I 1 / dE * dE l d( * \ M *\ 2 A s N) (/"> /»> /*' U) <V (A-35) 

fj, 

The phase space factor reads 

Af } (/at, //, ft, f q ) = - (1 + e ^ )(1+ \t t f (e t^ / + e ^)' (A.36) 

where we have used energy conservation, and Fermi-Dirac statistics for the leptons and quarks. 

The integration over q can now be performed analytically, reducing the dimensions of the 
collision integrals to two. These final integrals must be evaluated numerically and then summed to 
give CsMnI 

Cs, s [In] = CP + + CP + C<*> + C< B > + C< 8 > . (A.37) 
The integrals cjg' s '"' 6 ' ) are as follows: 



The first integral comes from evaluating the function Q\ a , and we have defined q = q/T: 

4 .» = a^i r * f +£ " n Af » ,« (A . 38) 

2 b ir s £ N y N v 2 J yN J 
r £l+m 7 A£n+£i) 2 -z 2 -q 2 



4 1} =/ dq- 



Isi-vn ' (£ N + £i) 2 -q 2 
Ay N (8i + £ N ) + z 2 log 

2 (Si + S N ) 



(£n—Vn) (2g;+£jy-yjy) 
(£n+Vn) (2£i+£ N +y N ) 



(A.39) 
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• The second integral comes from evaluation of O^: 



C 



— ^ ; — , (IS, / dS t lf\ 

s ' s 2^£ N y N v 2 J 7o 



(2) 



3T hfMrfu 



Jvn-£i 



2 2-2 
2T - <T 



4f, (£, + £jv) + log 



2 (£/ + £jv) 



(A.40) 



(A.41) 



The third integral comes from the Q,2a,i term: 



C, 



(3) 



3 T hfMnn 



rh (£n+vn) , , , , 



s ' s 2 6 vr 3 £: i vyjv v 2 

Jsi~y N (£ N +£i) 2 -q 2 

2 (Si + S N ) (£ N -2£ t + y N ) + z 2 log 



£t{2£i+£ N -y N ) 
(£ l +£ N -£t) (£n+Vn) 



2 (Si + £ N ) 



(A.42) 



(A.43) 



Integral four originates from the f^a a term: 



C 



(4) 



3T h 2 Mfhi 



s ' s 2 6 7r 3 £ N y N v 

1^ = - 



d£i I dStAWij.*), (A.44) 

J\ (2£i+£ N -y N ) 



2 ^-^ (£ N+£l ) 2 - z 2 -q 2 

dq g 

l£i-y N (£ N + £i) -q 2 

2 (Si + £ N ) (2£ l +£ N -2£ t - y N ) - z 2 log 



(£ l +£ N -£t)(2£ l +£ N -y N ) 



£t (£n+vn) 



2 (Si + S N ) 



(A.45) 



• Integral five relates to fl2b,i- 

(5 ) _ 3T h 2 t Mih x rvN 



C 



s,s 267r 3 £ 
4 5) = - / 

J V 



NVN V z 
dq ■ 



r\ (2£i+£n-Vn) , % , s 

o Jo 



(A.46) 



y N -£i (£n + £i) -q 2 

2 (Si + £ N ) (2£ t +£ N -2£ t - y N ) - z 2 log 

2 (Si + Sn) 



(£i +£ N -£ t ){2£i +£n-Vn ) " 

£t(£N+VN) 



(A.47) 
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Finally, the sixth integral derives from Cl2b,u- 



' 2 6 7T 3 ^ V 2 7 -A (£„+„„) 



2 (Zn+Vn) 
2£ '- £ > (£ N + St) 2 - z 2 - q 2 



4 6) = - / dg 

Jvn-£i 



(£ N + £ x f - q 2 
2 (£, + 5jv) (fjv - 2£ t + yjv ) - z 2 log 



(£i+£ N ~£t) (£n+Vn) 
£t (2£ t +£ N ~y N ) 



(A.49) 



2 (£j + f jv) 
A.1.2 Lepton asymmetry 

The s-channel collision term in Eq. (4.22) for tracking the lepton asymmetry is 

C sMa = ^j II (^f^ 2 ^^ 

(A.50) 

with phase space factor 

/at, /t, *,) = /,_! (/at (A + /, - 1) - /*/,) , (A.51) 

and matrix element given by Eq. (A. 11). 

Again we choose an explicit coordinate system and use rotational invariance 

q= |q| (o,o,i), 

pi = Ei (0, sin rj, cos r/) , (A.52) 
Pt = E t {coscj) sin i?, sin </> sin??, cost?). 

Following the same procedure as in section A. 1.1, we reduce the collision integral (A.50) to 

Cs,s[fi-j] = 27?r 3 E 2 J d cos tfd cos r]dE t dE N dq\M s \ 2 (A.53) 
/ Ef-E% + M+\q\* \J E 2 -(E l + E N -E t ) 2 + H 2 \ 

X I COS 7] — — — COS 17- 



2|q|£?i / V 2\q\E t J 

x a?-*) (/,_ 7 , f N , f t , t q ) e (^) e (E t ) e (s, + e n - E t ) , 

with phase space element 

using as usual energy conservation, and Fermi-Dirac statistics for the leptons and quarks. 

In analogy to section A. 1.1, we further reduce the collision integral (A.53) to a sum of six 
integrals with distinct integration ranges, 

CsAU = Cgj + + Cgj + C7g + Cgj + ^ (A.55) 
to be integrated numerically over two remaining degrees of freedom. These integrals are: 
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• First integral 

^ jf^ *» sr *> ^ * 

where is given by Eq. (A.39). 

(2) 

• Second integral with I s given by Eq. (A.41): 

(3) 

• Third integral (Is given by Eq. (A.43)): 



• Fourth integral (1^ given by Eq. (A.45)): 



(5) 

• Fifth integral (7^ given by Eq. (A.47)): 



• sixth integral (1^ given by Eq. (A.49)): 



A.2 t-channel 

A.2.1 Right-handed neutrino 

The collision integral for the t-channel process appearing in Eq. (4.20) is given by 

C S ,t[fN] = ^ J II jM^{^)^\PN+V q -Vt-Pl)\M t \ 2 k\ N) {f N J q J U ft), 
i=l,q,t 

(A. 62) 

with phase space factor 

A t (A ° (fN,f q Ji,ft) = [(1 - /iv)(l - f q )ftfi - f N f q {l - ft)(l - h)] , (A.63) 

and 



|M| 2 = 24/^™» (A.64) 
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is the matrix element. 

The reduction of the collision integral proceeds in the same way as for the analogous s-channel 
collision integral in section A.l. We use the momentum 

k = pi - pat = p q - pt, (A.65) 

and coordinates 

k = |k| (0,0,1), 

p q = E q (0, sin rj, cos 77), (A.66) 
PN = |Piv|(cos (j> sin#, sin <f> sini9, cos 1?). 
The matrix element in these coordinates reads 

» 2 (En - Eif - |k| 2 

Averaging over the incoming RHN direction and integrating over all azimuthal angles leads to the 
following 5-dimensional integral: 



CsA/n] = o7 o [ d cos fid cost] dE q dEid\]i\ \M 

E N p N J 



2 

t\ 



J , {EN + Eg-Etf-El-^ f E?-E 2 N + M 

x 8 cos 77 H <5 cos ■& - 



2\k\E q J V 2|k||pAr| 



x A ( t N) (f N , f q , fi, f t ) 6 (Eq) 6 (E t ) G (Ejv + E q - , 



(A.68) 



with 



.(«),, „ , e £ '+ £ * (-l + /jv + e £jv /iv) 

(1 + e^-ij (1 + e £ 9j (e^ + e £!+t< ?) 
assuming thermal equilibrium for the standard model particles. 

The integrals over cos 77 and cos 1? in Eq. (A.68) can be readily performed in the same manner 
as before, and, in the process, integration limits are derived for the integral over k = |k|. As 
it turns out, there are two possible lower limits for the ^-integration: fc m i nj i = E^ — E\ and 
&min,2 = E\ — p^. Here, a comment on the infrared cut-off is in order. In the integrated picture, 
a divergence occurs in the integral over t at t = 0, whose regulation requires the introduction of a 
Higgs mass in the propagator, i.e., \Mt\ 2 oc 1/t — > \Mt\ 2 oc l/(t — m\). In the full treatment, 
the matrix element has the form \Mt\ 2 oc 1/((En — Ei) 2 — \k 2 ), so that the equivalent divergence 
occurs in the integration over k at k = E^ — E\, i.e., at k = /c m in,i- This divergence can be avoided 
simply by modifying by hand the integration limit f k i — > J k 1+m<s> - There are no changes for 
those integrals with k m i n = /c m in,2- 

It is also possible to regulate the divergence by introducing a Higgs mass in the propagator, 
such as in the integrated picture. This modifies the integration over k not only for En > E\ 
(i.e., £; m i n = fc m i n) i), but also for E\ > -p^ (i.e., k m - m = fc m in,2)- In physical terms this procedure 
corresponds to giving the Higgs particle a mass whose magnitude can vary from zero up to possible 
thermal contributions, i.e., < m$/M < OAT/M, where m$(T) ~ 0.4 T is the thermal Higgs 
mass [6]. In the temperature regime relevant to leptogenesis, electroweak symmetry is unbroken 
and therefore leptons are massless. 
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Since we have so far not included thermal corrections in the present work, for consistency we 
prefer not to use the thermal Higgs mass. Furthermore, in a full thermal treatment, RHN decay into 
a lepton and Higgs pairs become kinematically forbidden at high enough temperatures, and the 
decay of a Higgs particle into a neutrino and lepton pair becomes viable [6]. Thus, in addition to 
determining the value of the infrared cut-off there is also the question of its interpretation. In view 
of these issues, we choose to deal with the infrared divergence using the simpler method of cutting 
off the integration over k at k m \ n = fc m in,i + m$, with ah = m$/M = 10~ 5 . 

After integrating over cos 77 and cos the original integral (A.68) is now split into four parts 



Cs,t[f~ 



(A.70) 



where the constituent integrals are as follows: 



First integral (with k = k/T): 



C 



(i) 



3T h 2 Mm x fi( £ N+y N ) 



S,t 



r(l) 



2 6 7T 3 SnVn v 2 Ji ( £ N _ yN+ahz) 

2£ q+ e N -e l {£n _ £) 2 _ z 2_~ e 
dk ■ 



dS, 



U2£ l ~£ N +y N ) 



r\{2£ 
J \a h z 



6Z q ^l?\ (A.71) 



/; 



-£i+a h z (£n ~ £l) ~ k 2 



2 (£ N - Si) (2£ q - a h z) - z 2 log 



(£ N +£ q -£i) a h z 



£ q (2(£ N -£ l )+a h z) 



2 {£ 



N 



• Second integral: 

y(2 ) _ 3T h 2 Mm x fh( £ N+y N ) 



C 



s,t 



r(2) 



2 6 7r 3 £ N y N 



d£i 

\{£ N -y N +a h z) - , 



poo 



d£ q A^4 2 \ (A.72) 



£ ' +VN dk {£N ~ £l)2 ~ z2 ~~ k2 



£ N -£i+a h z (£n ~ £l) ~ k 2 

2 (£ N - £1) (2£i - £ N + y N - a h z) - z 2 {log 



■£n+Vn) 



(£ N -2£i-y N ) {2{£ N ~£ l )+a h z) 



2 (£ N - £1) 



• Third integral 

,(3) :'•/ hf^l'in r " 



C 



s > f 2 6 vr 3 £ N y N v 2 



POO f 

/ d£i \ 

JWn+Vn) 

z 2 



r(3) 



2£ q+ £ N -£ l {£n _ £i) 2 _ z 2_~ k 2 

Ufv 



\(2£ l ~£ N +y N ) 



(2£ l ~£ N -y N ) 



d£ A (N) I (3) 

ac q J H 1 t 1 



£i-y N (£ N -£i) -k 2 

2 (£ N - Si) (£ N + y N + 2(£ q - £ t )) + z 2 (log 

2 (£ N - £1) 



(A.73) 



-£g(£N~yN) 

(£ N +£ q ~£ l ) (£ N ~2£ l +y N ) 
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• Fourth integral: 

, 4)= 3T jjMg,r r» j.^^ (A . 74) 



£i-yjv (E N -Ei) -k 2 

4 (£ jy - £/) + z 2 log 



(£n-Vn) (£n-Vn-2£i) 

(£n+vn) {£n+vn-2£i) 



2 (£tv - £i) 

A.2.2 Lepton asymmetry 

The i-channel collision integral for the lepton asymmetry evolution is 

(A.75) 

with 

Af _I) /*, /,) = h-j (f q (ft + f N -l)~ ftf N ) ■ (A.76) 
The matrix element is the same as for the RHN given in Eq. (A.64), and we use the momentum 

k = pat - p/ = p q - p t , (A.77) 

and coordinates 

k = |k| (0,0,1), (A.78) 
p q = E q (0, sin 77, cos r/), 
pi = Ei (cos (j) sin 1?, sin 4> sin 1?, cos #) . 

Following the method of the previous sections, we reduce the integral (A.75) to 

Cs,t[fi_j] = 2?7r 3 E 2 J d cos fid cos T]dE q dE N d\k\ \M t \ 2 

J + E q - E N ) 2 - E 2 - \k 2 \ \ J E 2 N -E?-M 2 -\k 2 \ 
x 6 cos rj H <5 cos $ 



2|k| ^ y v 2|k|£, 

(A.79) 



x a? I} / t , /at, /,) e (e ? ) e (E^) e + E q - e n ) , 



with 

Ai-D (f ffAif] _ f e£q (e £l i-^ + fN)-e^f N ) 
* /*' /*> " /i-i (i + e f 9 ) (gfiv + e^+e,) (A - 80) 



as the phase space factor. 

Integrating over cos 77 and cos we split up the integral (A.79) into four parts, 

Cs,t[fi-i\ = C { s} + 4 2 J + Cs} + (A.8D 

with: 
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• First integral: 

= /-^^^ (fly - gQ 2 - Z 2 - fc 2 

Je N -e,+a h z (S N - Si) 2 -k 2 



-£i+a h z (£n - Si) - k 2 
2 (£ N - £^ (2 (£ N -£ q - £^ + a h z) + z 2 log 



-£n+£i) (2(£jv-^)+ a /i 2; ) 



2 («fjv - Si) 

• Second integral: 

4i = 6 3 3 j ^ r 2 . ^ r ^ Af- 7) i?\ (a.83) 



J£ N -£l+a h z {£n - £i) - k 2 

2 (.Sat - £j) (2£; -£ N + y N - a h z) - z 2 log 



{£n+Vn) ah* 



{2£i~£ N +y N ) (2{£ N -£ l )+a h z) 



2 (gjv - £i) 

• Third integral: 



2 7r £nVn v 2 J z J\{£ n -vn) 

7 (3) = ^ (^-^ Z 2_fc2 



2 - fj) (2£ g - £tv + yjv) + ^ 2 log 



{£ N —£ q —£l) (£n—Dn) 



£ q (£n~2£i+vn) 



2 (£ N - Si) 

Fourth integral: 



|(f]V+J/Jv) 



(Sn-Si) - k 2 



I£i-vn {Sn-Si) - k 2 

4 (.Sat - fj) yAr + z 2 log 



(£ N ~2£i~y N ) (Sn-Un) 
(£ N -2£ l +y N ) (Sn+Vn) 



2 (gjv - 

B. Evolution of the top Yukawa coupling 

To determine the gauge and Yukawa couplings at some energy scale fi = log (T/mz), we use the 
renormalisation group equation, 

^"8^- (B - 1} 
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where % denotes the corresponding gauge group of the standard model. The constants for the gauge 
couplings are 



. /41 5 16 _ 

(<*><*><*) = (j0- 3' "9 " 7 

At one loop the solution of (B.l) for the gauge couplings yields 



(B.2) 



9i 
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Neglecting contributions from the bottom and charm Yukawa couplings, the renormalisation group 
equation for the top Yukawa coupling at one loop is given by [34] 



dfi 8vr 2 
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A detailed study of the evolution of quantities relevant for leptogenesis can be found in [35]. 
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